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METHOD AND ARRANGEMENT FOR NOISE VARIANCE AND SIR 

ESTIMATION 

5 Field of the Invention 

This invention relates to noise variance and 
signal /interference ratio (SIR) estimation, and 
particularly though not exclusively to such estimation in 
10 wireless communication receivers. It will be understood 
that, as used herein, the terms 'noise' and 
'interference' are to be considered synonymous, with each 
encompassing both noise and interference. 

15 

Background o£ the Invention 

In the field of this invention it is known that many 
parts of a wireless communications receiver often require 
20 an estimation of noise variance and/or SIR. This is 
needed for purposes of power control, threshold 
determination for various algorithms, quantisation of 
soft -decision information for channel decoding purposes 
to name but a few. 

25 

For BPSK (Binary Phase Shift Key) and QPSK (Quadrature 
Phase Shift Key) modulation the conventional method for 
estimating the SIR at the output of a detector relies on 
estimating output noise variance using the following 
30 equality known (for example) from the publication by 
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Papoulis and Pillai, entitled ^Probability, Random 
Variables and Stochastic Processes', 3rd Ed. 1991, 

a,'=£(|rf<'f)-4?W|)' 

where represents variance, E represents mean value 
5 and "rf^*^ represents the detector output J 

This yields the following result: 

where SIR represents the SIR of the sequence at the 
10 detector output, and represents the average power of 
the sec[uence at the detector output. 

However, this approach has the disadvantage (s) that the 
accuracy of this method at low SIR is poor since it 

15 suffers from a bias term. An analysis of the bias term 
and a correction method has been suggested in GB patent 
application no. 012 8475.1. However, the suggested 
correction method requires a look-up table to correct for 
the aforementioned problem, and the estimation variance 

20 is also increased when correcting the bias. 

A need therefore exists for a method and arrangement for 
noise variance and SIR estimation wherein the 
abovementioned disadvantage (s) may be alleviated. 

25 
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Statement o£ Invention 

In accordance with a first aspect of the present 
invention there is provided a method for noise variance 
5 estimation as claimed in claim 1, 

In accordance with a second aspect of the present 
invention there is provided an arrangement for noise 
variance estimation as claimed in claim 10, 

10 

Preferably, the second noise variance signal is produced 
by applying to the first noise variance signal a function 
substantially equal to the detector's transfer function. 

16 Preferably, the first noise variance signal is derived 
from a midamble portion of the received signal . 

Preferably, an estimate of total power at the detector 
output is produced from the second noise variance signal 
20 and an SIR signal representative of SIR in the received 
signal . 



Brief Description of the Drawings 

25 

One method and arrangement for noise variance and SIR 
estimation incorporating the present invention will now 
be described, by way of example only, with reference to 
the accompanying drawings, in which: 
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FIG. 1 shows a block schematic diagram illustrating 
a 3GPP radio coiranunication system in which the 
present invention may be used; 

5 FIG. 2 shows a block schematic diagram illustrating 

. a generic burst structure employed in the system- of 
FIG. 1; and 

FIG. 3 shows a block schematic diagram illustrating 
10 a detector arrangement incorporating a noise 

variance and SIR estimator utilising the present 
invention. 



15 Description of Preferred Eznbodiment ( s ) 

The following preferred embodiment of the present 
invention will be described in the context of a 3GPP (3^^ 
Generation Partnership Project) UMTS (Universal Mobile 

20 Telecommunication System) Radio Access Network (UTRAN) 
system operating in TDD mode. Referring firstly to 
FIG. 1, a typical, standard UTRAN system 100 is 
conveniently considered as comprising: a terminal/user 
equipment domain 110; a UMTS Terrestrial Radio Access 

25 Network domain 12 0; and a Core Network domain 130. 

In the terminal/user equipment domain 110, terminal 
equipment (TE) 112 is connected to mobile equipment (ME) 
114 via the wired or wireless R interface. The ME 114 is 
30 also connected to a user service identity module (USIM) 
116; the ME 114 and the USIM 116 together are considered 
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as a user equipment (UE) 118, The UE 118 coinmunicates 
data with a Node B (base station) 122 in the radio access 
network domain 120 via the wireless Uu interface. Within 
the radio access network domain 120, the Node B 122 
5 communicates with a radio network controller (RNC) 124 

via the lub interface-. The RNC 124 communicates with other - 
RNC's (not shown) via the lur interface. The Node B 122 
and the RNC 124 together form the UTRAN 126. The RNC 124 
commxinicates with a serving GPRS service node (SGSN) 132 

10 in the core network domain 130 via the lu interface. 
Within the core network domain 130, the SGSN 132 
communicates with a gateway GPRS support node (GGSN) 134 
via the Gn interface; the SGSN 132 and the GGSN 134 
communicate with a home location register (HLR) server 

15 136 via the Gr interface and the Gc interface 

respectively. The GGSN 134 communicates with public data 
network 138 via the Gi interface. 

Thus, the elements RNC 124, SGSN 132 and GGSN 134 are 
20 conventionally provided as discrete and separate units 
(on their own respective software/hardware platforms) 
divided across the radio access network domain 120 and 
the core network domain 130, as shown the FIG. 2. 

25 The RNC' 124 is the UTRAN element responsible for the 

control and allocation of resources for numerous Node B's 
122; typically 50 to 100 Node B's may be controlled by 
one RNC. The RNC also provides reliable delivery of user 
traffic over the air interfaces. RNC's communicate with 

30 each other via the lur interface. 



BEST AVAIUBLE COPY 



wo 2005/015790 



6 



PCT/GB2004/003368 



The SGSN 132 is the UMTS Core Network element responsible 
for Session Control and interface to the HLR. The SGSN 
keeps track of the location of an individual UE and 
performs security functions and access control. The SGSN 
5 is a large centralised controller for many RNCs. 

The GGSN 134 is the UMTS Core Network element responsible 
for concentrating and tunnelling user data within the 
core packet network to the ultimate destination (e.g., 
10 internet service provider - ISP) . 

Such a UTRAN system and its operation are described more 
fully in the 3GPP technical specification documents 
3GPP TS 25.401, 3GPP TS 23.060, and related documents, 
15 available from the 3GPP website at www.3gpp.org, and need 
not be described in more detail herein . 

The physical layer of UTRA TDD mode provides physical 
channels that carry transport channels from the MAC 

20 (Medium Access Control) sub-layer of UMTS Layer 2. A 
physical channel is defined by frequency, timeslot, 
channelisation code, burst type, and radio frame 
allocation. In UMTS Layer 2, in each time slot, three 
burst structures (as shown generically in FIG. 2) can be 

25 supported, each consisting of two data fields (210, 230), 
a midamble (220) and a guard period (240) . 

The data fields contain the data symbols from the 
transport channels, after the processes of coding, 
30 multiplexing, interleaving, and modulation. The midamble 
field contains the training sequence, which is used in a 
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number of Layer 1 algorithms, such as channel estimation . 
The guard period, GP, is used to accommodate any timing 
inaccuracies, from propagation delays, channel 
dispersion, and power ramping in the transmitter. The 
5 different burst types and their associated field lengths 
in -chips are given in the table below: 



Burst Type 


Data 
Field 1 


Data 
Field 2 


Midamble 

m 


GP 


Burst Type 1 


976 


976 


512 


96 


Burst Type 2 


1104 


1104 


256 


96 


Burst Type 3 


976 


880 


512 


192 



The received sequence in the data payload areas of the 
burst is given by 

e= Ad+n 

where 

d = (c,0)r,d(^)^ fS^'^J ={d,.d,....,d^y . 

(.)'' denotes transposition, K is the number of data 
sequences k-\,...,K, and N is the number of symbols per 
data sequence n-\,...,N . The data sequence for the A:** user 

is given by d(*^ = (jf*\</?^---»'='wT • '^^^ noise sequence n, 
denoted by n = (n„nj,...,/i;^+B,_,)r , has zero mean and covariance 
matrix R„=£(nn"), where denotes conjugate 
transposition. The matrix has dimensions {NQ + W -\)y.KN 
and the elements are given by: 

_rb^*^ A = 1...A:, l = \...Q-\-W-\,n = \...N 
0 else 



15 
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where b^*^ =h^*^*C^*^, h^*^ is the impulse response of the 
user, W is the length of the channel impulse response, * 
denotes discrete time convolution, C^*^ is the spreading 
code of the user and Q is the length of the spreading 
5 sequence in chips. 

The output of the detector is given by 

d = /(e)-r+z 

where /(.) denotes the transfer function of the detector, 
10 the vector r contains the desired symbols, and the vector 
Z contains noise plus interference. The average power for 

the sequence at the output of the detector is given by 

15 Expanding P^*^ produces 

Under the assumption that the noise is uncorrelated with 
the signal vector r, the average power for the k^ 
20 sequence becomes 

where E{) is the statistical average, cr^^ =^|^„+(*-i)0f ) is the 

noise variance at the output of the detector, and 

^|';+(*-i)i2f ) ^® signal power for the k*^ sequence. The 

25 SIR at the output of the detector for the k^ sequence is 
therefore given by 
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The conventional method for estimating the SIR relies on 
estimating the detector output noise variance using the 
5 following equality mentioned above: 



10 As discussed above, the accuracy of this approach at low 
SIR is poor since it suffers from a bias term, which may 
be corrected by use of a look-up table. 

As will be discussed in greater .detail below, the 



do not suffer from such a bias term and therefore do not 
require a look-up table to correct for the aforementioned 
problem. 

20 Referring now to FIG. 3, in detector arrangement 300 

(which may be used in Node B 122 or User Equipment 118) a 
new technique for estimating the SIR at the output of a 
detector 310 uses an estimate of the detector input noise 
variance (derived in channel estimator 320) to provide an 

25 estimate of the detector output noise variance (derived 
in SIR estimator 330) . In the channel estimator 320 the 
detector input noise variance is derived from the 




to yield the following result: 




15 following preferred embodiments of the present invention 
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midatnble portion of the burst using locally reconstructed 
replicas of the transmitted signal which are removed from 
the received set of samples providing a residual term 
that contains the noise of the channel . In the SIR 
5 estimator 330, by deriving the transfer function of the 

detector an -estimate of the detector output noise- - 

variance is estimated. 

The estimated output noise variance then allows an 
10 improved estimate of the SIR ( SIR^^K.SIR^^b at the detected 
output. Typically, the SIR at the output of the detector 
is used for deriving soft decision cjuantisation levels 
for application to channel decoding algorithms. 

15 In the following description, two types of CDMA (Code 

Division Multiple Access) detector are considered, namely 
single user detector (SUD) and multiuser detector (MUD) • 
It will be understood that the invention is also 
applicable to other types of detector such as a RAKE 

20 receiver . 

The technique described here is based on first estimating 
the noise variance at the input to the detector and then 
mapping the input noise variance to the output noise 
25 variance using the transfer function of the detector. 

The process of estimating the noise variance at the input 
to the detector is carried out using the midamble portion 
of the burst- Considering the received sequence of chip 

30 spaced samples e = (e,,e2,-.-,e^^)^ / where is the burst 



^EST AVAILABLE COPY 



wo 2005/015790 PCT/GB2004/003368 

11 



length, a locally reconstructed version for the midamble 
portion of the burst is defined as follows: 

5 = (^.^ 5^.,J=(E('>'*''>"'"j]'' 

where denotes locally reconstructed version and m^*^ is 
5 the midamble sequence of the A"* user. The noise variance 
estimate at the input to the detector is given by 

where X<L^ and the starting position is W since the 
first W-\ samples from the midamble portion of the burst 
10 are corrupted by the data portion of the burst. 

Multiuser Detection 

Under the assumption that the noise is white with 
15 variance cr^ , the Minimum Mean Squared Error (MMSE) block 
linear equalizer solution to symbol estimation is given 
(as known from the publication of Klein, Kaleh and Baier 
entitled *Zero Forcing arid Minimum Mean -Square -Error 
Equalization for Mutliuser Detection in Code-Division 
20 Multiple-Access Channels' in IEEE Trans VT, VOL.45, No. 2, 
May 1996, pp276-287) by 

a = (a" A + A^e = /(e) = r + z 
where 1 is the identity matrix and 

r = (A^A + aM)r'A^Ad = (r„r„...r^f (D 

26 z = (A^A+cT^l)"^A''n = (z„z,,...,z^J^ (2) 
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From equation (2), the noise variance seen at the output 
of the MUD is given by 

where H denotes vector norm, and cr^ represents the noise 
5 variance at the input of the MUD. By replacing with 
the estimate of the MUD input noise variance a', we have 
a direct method for estimating the MUD output noise 
variance &l . For completeness, the estimate &l of the MUD 
output noise variance can be written as 

where f(&^) represents the noise transfer function of the 
detector. 

Using the new estimate for the output noise variance, the 
15 SIR at the output of the MUD for the k'^ sequence is 
defined by 

where the error term is given b^5tb I AVAILAbLt COPY 

2 ^2 

20 

It is clear from the above set of equations that when 
-<y^ we have the following 
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(*-l)g| ) 



It will therefore be understood that the accuracy of the 
above technique is directly related to the quality of the 
5 noise variance estimate, , at the input of the MUD. 

Single User Detection 

For the single user detector case the received sequence 
10 is written as 

e=Ad+n=HCd+n 



The matrix H has dimensions {NQ-hW-ljxNQ and its 
elements are given by 

[ 0 else 

where h = {h„h„...,Kj , i = l,...,NQ + W -1 , and v = l NQ . The 

matrix C has dimensions NQxKN and its elements are 
given by 



15 



<^e(i.-i)f9.Ar(*-iK» I Q g/^g 

20 

For Minimum Mean Squared Error (MMSB) symbol estimation 
and under the assumption that the noise is white with 
variance <t^ , the output of the SUD is given by 

d = C" (h"H + a^y H"e = /(e) = r + z 

25 where 

r = C"(H"H+o-^l)r'H*HCd = (/'„/2,.../V;,)'' 
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z = (h^H + cr ^l)r' H'^n = (z„ z„ . . . , z^^ / 



(4) 



From equation (4) , the noise variance seen at the output 
of the SUD is given by 

- - ■ <T.'=4.<.-.tor)=ex|hr-.<.'y'ihiiv= 

where ||| denotes vector noinna, the multiplier G comes from 
the matrix C, and in general G = ||c^*^f =0. 



By replacing cr^ with the estimate of the SUD input noise 
variance , we have a direct method for estimating the 
SUD output noise variance a] . For completeness, the 
estimate &] of the SUD output noise variance can be 
written as 

15 where G is replaced with Q and /(o'^) is the noise 

transfer function of the detector. Using the new estimate 
for the output noise variance, the SIR at the output of 
the SUD for the sequence is defined by 



20 where the error term s[&^) is given by 

*(<T=)=cxihr 



2 * 2 



It is clear from the above set of equations that when 
-a^ we have the following 
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It will therefore be understood that the accuracy of the 
above technique is directly related to the quality of the 
5 noise variance estimate, , at the input of the SUD. 

It will be appreciated that the method described above 
for noise variance and SNIR estimation may be carried out 
in software running on a processor (not shown - e.g., in 
10 User Equipment such as 118 or a Node B such as 122) , and 
that the software may be provided as a computer program 
element carried on any suitable data carrier (also not 
shown) such as a magnetic or optical computer disc. 

15 It will be also be appreciated that the arrangement 

described above for noise variance and SNIR estimation 
may alternatively be carried out in hardware, for example 
in the form of an integrated circuit (not shown) such as 
an PPGA (Field Programmable Gate Array) or ASIC 

20 (Application Specific Integrated Integrated Circuit) . 

It will be understood that the method and arrangement for 
noise variance and SIR estimation described above 
provides the following advantages that the accuracy of 
25 this technique is not poor at low SIR, since it does not 
suffer from a bias term, nor does it require correction 
therefor using a look-up table. An additional advantage 
is that any increase in estimation variance resulting 
from bias correction may be avoided. 

30 
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